Charmonium spectroscopy with heavy Kogut-Susskind quarks 
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Charmonium spectroscopy with Kogut-Susskind valence quarks are carried out for quenched QCD at f3 = 6.0 
and for two-flavor full QCD at j3 = 5.7. Results for IP— IS mass splitting and estimates of a—!—(mz) are reported. 
Problems associated with flavor breaking effects and finite size effects of IP states are discussed. 



1. INTRODUCTION 

An interesting application of lattice heavy 
quark QCD is a determination of the strong cou- 
pling constant through the IP-IS mass splitting. 
A shortcoming of the initial studies[l], however, 
was the use of quenched simulations, necessitat- 
ing an estimate of sea quark effects via the per- 
turbative renormalization group. While a mean- 
field estimate of the renormalized coupling in 
quenched and full QCD provided support for the 
validity of this procedure [2], a direct quarkonium 
spectrum simulation in full QCD is clearly desir- 
able. An exploratory calculation has been carried 
out in Ref.[2] using Wilson valence quarks on full 
QCD configurations generated with two flavors 
of Kogut-Susskind (KS) quarks. Here we report 
on our attempt to employ the KS action also for 
heavy valence quarks. Our study is carried out in 
two steps; a detailed quenched study of the char- 
monium spectrum for examining the viability of 
the KS action for heavy quarks in relation to fla- 
vor breaking effects, followed by an extension to 
full QCD and estimates of the strong coupling 
constant. 
* presented by S. Aoki 



Table 1 

Run parameters 



/3 



(GeV) 



#conf. 



m h 



quenched 6.0 1.88(6) 24 3 x 40 



16- 5 x 40 



full QCD 5.7 2.23(9) 20 3 x 20 
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2. SIMULATION 

Parameters of our runs are listed in Table 1 in- 
cluding the values of heavy quark mass m/ja and 
the lattice spacing determined from m p . For full 
QCD we use the configurations generated previ- 
ously with the sea quark mass of m sea a = 0.01 
and 0.02[3], taken at 20 trajectory intervals and 
periodically doubled in time. Quenched configu- 
rations are those employed in Refs. [4,5]. 

We employ 128 meson operators extended over 
a 2 4 hypercube and non-local in time. These are 
classified into 8 GTS irreducible representations 
for 7) c , 12 for J/i}>, 6 for h c , and 4 and 6 for Xc 
and Xcj6]- Eight wall sources at t = and 1 are 
used to maximize propagator signals [4] calculated 
with Landau gauge fixing. Hadron propagators 
are averaged over members of irreducible repre- 
sentations before fitting. Errors are estimated by 
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Figure 1. Operator dependence of quenched me- 
son masses at (3=6.0 and m/ja=0.3 on a 24 3 X 40 
lattice. 

the single elimination jackknife procedure. 

3. CHARMONIUM SPECTRUM 

The operator dependence of charmonium 
masses is illustrated in Fig. 1. A notable feature 
is that 7j c (NG) in the Nambu-Goldstone chan- 
nel is much lighter than the others (rj c (non-NG)) 
even at heavy quark masses. It is quite disturbing 
to find that the mass splitting increases for larger 
quark mass, as shown in Fig. 2 where the splitting 
relative to the Nambu-Goldstone pseudo scalar is 
plotted as a function of m^a for the 5 GTS ir- 
reducible representations for non-NG r\ c . Other 
states also exhibit a similar increase of mass split- 
ting with m^a, though much smaller in magni- 
tude. 

There is no reason to favor a particular channel 
when comparing measured masses with the ex- 
periment. We, therefore, consider meson masses 
averaged over GTS representations having the 
same continuum J pc weighted by their dimen- 
sions. For j; c , however, we separate out rj c (NG). 
In Fig. 3 we plot the averaged masses in physical 
units for quenched QCD employing the scale de- 
termined by m^, together with the experimental 
values. For the states other than r\ c the calcula- 
tion reproduces qualitative features of the char- 
monium spectrum. However, we find mj/j, — 
m ri c (NG) = 404(4)MeV, which is much larger than 
the experimental value 118MeV, while a factor 
two smaller value of 51(6)MeV is obtained with 
the average over non-NG r\ c . 

Results for the charmonium spectrum for full 
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Figure 2. m Vc a — m v ^ NG ^a for 5 non-Nambu- 
Goldstone r\ c states as a function of m/j<x at /3=6.0 
on a 24 3 x 40 lattice in quenched QCD. 

QCD is quantitatively similar to the quenched re- 
sults. For J ftp — 7] c mass splitting, in particular, 
we find 442(4)MeV for rj c (NG) and 62(7)MeV for 
7j c (non-NG). A 15% difference in the lattice spac- 
ing between our full and quenched simulations 
makes it difficult to see if the larger splitting for 
full QCD represents a physical effect of sea quarks 
that is expected from renormalization group con- 
siderations. 

4. IP-IS MASS SPLITTING 

Our results for the charmonium IP-IS mass 
splitting for quenched and full QCD are summa- 
rized in Fig. 4. For quenched QCD the splitting is 
almost independent of m^a for each lattice size, 
which is consistent with the experiment. There 
exists, however, a 20% difference between the val- 
ues obtained with rj c (NG) and 7] c (non-NG). 

Another important point is a significant lattice 
size dependence of the splitting between L = 16 
and 24 observed for quenched QCD. We found 
that this arises from a size dependent mass of the 
IP state h c , which increases by 2.4(9)% between 
the two lattice sizes (for J/ip and rj c the change 
is less than 0.2 %). Though small for the mass 
itself, this shift translates into a 20 % increase for 
the IP— IS splitting. 

The physical lattice size of La ~ 1.8fm for 
full QCD on an L = 20 lattice is comparable 
to La ~1.7fm for quenched QCD on an L = 16 
lattice. In view of the size dependence for the 
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Figure 3. Averaged meson masses in physical 
units as functions of m/j<x at /3=6.0 on a 24 3 X 40 
lattice in quenched QCD with the scale fixed by 
vn p . Crosses correspond to the point where the 
mass of J ftp agrees with the experiment. 

quenched result it is likely that the full QCD 
value for the IP-IS mass splitting is an underes- 
timate. In fact, employing the experimental value 
Ato 1p _i S = 458MeV, we find a - 1 =2.70(20)GeV 
with 7] C (NG) and 3.59(36)GeV with rj c (non-NG) 
at m sea a=0.01 and ?7i^a=0.3, which are signif- 
icantly larger than the value 2.23(9)GeV esti- 
mated from m p . The corresponding values for 
the quenched case on an L = 24 lattice are 
1.87(16) GeV with r) c (NG) and 2.32(25)GeV with 
7] c (non-NG) at m/j<x=0.4, which are compared 
with 1.88(6)GeV from m„ 



o p . 



.(5) 



5. DETERMINATION OF a^L-(m z ) 

Our method for estimating the strong coupling 
constant follows that of Ref.[l]. The a\r cou- 
pling is used for renormalization group evolu- 
tion and fj, = 0.4 — 0.75 GeV is taken for the 
charmonium scale. For quenched QCD we cal- 
culate a-—(mz) for three values of mi,a and 

MS\ ^' . n 

the two choice of inverse lattice spacing corre- 
sponding to 7j c (NG) and 7j c (non-NG). Averaging 
the resulting six values of a--j—(mz), we obtain 

a ( IL(m z ) = 0.111(4) (4). The first error reflects 
the ambiguity of the charmonium matching scale 
and the second the choice of 7j c (NG) or 7j c (non- 
NG) for determining the scale. For full QCD a 
similar analysis yields a——(mz) = 0.119(3)(3). 
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Figure 4. IP-IS mass splittings in lattice units 
for quenched (open symbols) and full (filled sym- 
bols) QCD for various quark masses as a function 
of 1/ La with the scale fixed by m p . 

The true value must be smaller than this due to 
an over-estimate of the lattice scale as discussed 
above. 

6. Summary 

Charmonium spectroscopy with the KS action 
suffers from a sizable flavor breaking effect whose 
magnitude increases for heavier quarks, partic- 
ularly for 7j c , at a currently used lattice spac- 
ing. Systematic uncertainties in a——(mz) due to 
this source is difficult to be diminished. Whether 
modification of the action is possible for reducing 
the effect is an open problem at present. 
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